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Abstract 

The two-body t-matrix is calculated directly as function of two vector mo- 
menta for different Malfliet-Tjon type potentials. At a few hundred MeV 
projectile energy the total amplitude is quite a smooth function showing only 
a strong peak in forward direction. In contrast the corresponding partial wave 
contributions, whose number increases with increasing energy, become more 
and more oscillatory with increasing energy. The angular and momentum 
dependence of the full amplitude is studied and displayed on as well as off the 
energy shell as function of positive and negative energies. The behavior of the 
t-matrix in the vicinity of bound state poles and resonance poles in the second 
energy sheet is studied. It is found that the angular dependence of T exhibits 
a very characteristic behavior in the vicinity of those poles, which is given 
by the Legendre function corresponding to the quantum number either of the 
bound state or the resonance (or virtual) state. This behavior is illustrated 
with numerical examples. 
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I. INTRODUCTION 



At low energies in the MeV and the few tenth of MeV region very few angular 
momenta contribute to the nucleon-nucleon (NN) scattering process. Consequently a de- 
scription using angular momentum decomposition is an adequate tool for carrying out scat- 
tering calculations. However, at intermediate energies, i.e. energies of a few hundred MeV, 
and higher energies very many angular momenta contribute to the scattering amplitude. 
There those individual contributions to the scattering amplitude at a fixed high angular 
momentum oscillate strongly in angle, whereas the total amplitude is much smoother. This 
suggests the direct determination of T as function of the initial and final momentum vectors 
avoiding angular momentum decomposition totally. For NN scattering investigations of this 
kind have already been undertaken [0-^. 

The choice of momentum vectors as adequate variables is also suggested from the NN 
force. The dependence on momentum vectors in the case of the widely used one-boson- 
exchange force is for instance rather simple, whereas the partial wave representation of this 
force leads to complicated expressions This is already apparent in the most simple case 
of a scalar meson propagator, l/((q' — q)^ — /U^), which is the central ingredient to any 
NN force. In a partial wave decomposition this is represented in the form -^Qi{z), where 

z = (q'^ + q^ + /i^)/2g'g and Qi{z) is the Legendre function of the second kind. For large 
values of I the latter requires some care in order to be handled numerically correctly. 

Scattering of more than two particles requires two-body t-matrices off-the-energy-shell 
(off-shell for short notation) as dynamical input, which are apparently easier to handle if 
they enter the calculations as smooth functions instead of strongly varying partial wave com- 
ponents. Of course those remarks also apply to the treatment of scattering processes of more 
than two particles at intermediate energies, which themselves are also treated more econom- 
ically and transparently using momentum vectors instead of partial wave representations. 
Calculations of three and more particles use as input fully off-shell two-body t-matrices, 
whose properties as functions of arbitrary initial and final momenta and in general positive 
and negative energies should be well understood. Specifically at negative energies there may 
be bound state poles and in the second energy sheet there may be poles related to virtual 
states and resonances. 

Our aim in this article is to generate two-body t-matrices directly in a three-dimensional 
form and display their properties as function of the magnitudes of the off-shell momenta, the 
angle between the two momentum vectors and of the energy. We are not aware of a similar 
study in this generality in the literature. Usually only partial wave projected amplitudes 
are displayed and discussed [^. In Section II we describe our solution of the two-body 
Lippmann-Schwinger equation directly as function of the momentum vectors and illustrate 
the on- and off-shell properties of the resulting t-matrices obtained with simple Yukawa 
type two-nucleon potentials. In Section III we discuss the pole structure of the t-matrix as 
function of the energy and illustrate its angular and energy behavior at and around bound 
state poles. In Section IV we discuss and illustrate the behavior of the t-matrix for virtual 
and resonant states in the second energy sheet. We conclude in Section V. 
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II. THE ON- AND OFF-SHELL TWO-BODY T-MATRIX AT POSITIVE 

ENERGIES 



Two-body scattering is governed by the Lippmann-Schwinger equation 

T = V + VGoT, (2.1) 

where V is the two-body (e.g. two-nucleon) potential, Gq = {z — Hq)^^ the free two- 
body propagator and T the transition operator. In momentum space its matrix elements 
T{q',q,z) = (q'|T(2;)|q) obey the integral equation 

r(q', q, z) = V(q', q) + | rf^g'V(q', q")-^r(q", q, z). (2.2) 

m 

Here q are the relative momenta, m the common mass of the two particles and z an arbitrary 
energy. We use a nonrelativistic framework. In this article we restrict ourselves to two 
spinless particles and local potentials. Therefore, V{q', q) as well as T(q', q, z) are scalar 
functions: 

V{c{,q) = V{q',q,^-ci) (2.3) 

and 

T(q',q)=T(g',g,q'-q). (2.4) 

In Eq. ( |2.4] ) we dropped the parametric dependence on the energy z. This notation then 
leads to the explicit form of Eq. ( p^ ) 

T{q', q, x') = V{q\ q, x') + / dq"q"^ / dx" / d^"V{q\ q", y) -^T(g", g, x"), (2.5) 

Jo J-i Jo z — ^ 

m 

where x' = q' ■ q, x" = q" ■ q, and y = q" ■ q'. We can express y through x' and x" as 

y = x'x" + Vl - a;'Vl - x"^ cos^" (2.6) 
where the arbitrary azimuthal angle (f for q is chosen to be zero. If we define 

;>oo 

v(q' , q, x' , x) = / d(pV(q, g, x'x + Vl — x'^Vl — x'^ cosip), (2.7) 
Jo 

the integral equation Eq. ( |2.5| ) becomes 

1 1"°° 1 

T{q\q,x') = —v{q\q,x\l)+ dq"q"' dx"v{q\q\x\x") T{q\q,x") (2.8) 

ZTT Jo J-1 2—2— 

m 

This is a two-dimensional integral equation in the off-shell momenta q' {q") and the cosine 
of the 'scattering angle' x' {x"). 

In this Section we consider the solutions of the integral equation Eq. ( |2.8|) at 

positive energies, i.e. we choose z = E + ie = ^ + ie, corresponding to the incoming 
momentum qo. In order to obtain insight into the behavior of the t-matrix, we will consider 
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the on-shell element T{qo,qo,x,E), whose square is proportional to the differential cross 
section as well as the half-shell, T(g, go, x, E), and fully off-shell, T(g, q', x, E) t-matrix. 

We solve the two-dimensional integral equation typically using 24 or 32 q-points 
and 24 x-points. The Cauchy Singularity is separated into a principal value part and a 5- 
function part, and the principal value singularity is treated by subtraction. The integration 
interval for the q-integration is covered by mapping the Gauss-Legendre points u from the 
interval (0,1) via q = ctan(|M) to the interval (0,cxd). Typical values of c are 1000 MeV/c. 

A very stringent test for our numerics is the off-shell unitarity relation, which is a 
direct consequence of Eq. (|2.2|). In our two-dimensional form it reads 



7J- rl /•2n 

ImT{q', g, x') = --mqo J ^ dx" dip"T{q', go, y)7'*(g, go, x") (2.9) 

where y is given in Eq. ( p.6|) . With Eq. ( p.9|) we allowed for the most general case of the 
unitarity relation, where the energy z = E = is not related to the incoming momentum 
g, thus q ^ q' go- In our numerical tests Eq. ( pl9|) was fulfilled for arbitrary g', g, and x' 
values with an accuracy below 0.001% with the above quoted number of integration points. 
As main application we choose potentials of the Malfiiet-Tjon R] type, i.e. 



V{r) = Vr Va , (2.10) 



and consequently 



In the case of a Malfliet-Tjon type potential the yj-integration of Eq. ( |2.7| ) can be carried 
out analytically with the result 



v{q',q,x',x) = - 

71 



V, 



R 



(2.12) 



(g'2 _l_ g2 _ 2q'qx'x + /ij?)^ — Aq''^q'^{l — a;'^)(l — 

Va 

_|_ g2 _ 2q'qx'x + fiAY ~ 4:q''^q'^{l — x'^)(l — x^) 

The parameters used for V are given as V^^^ in Table I. Note that they are slightly different 
from the ones used in Ref. 0. 

As first numerical example we would like to demonstrate the connection of the 
angle depended on-shell amplitude T{qQ,qQ,x,E) and its representation in terms of partial 
wave amplitudes, 

°° 2/ + 1 

T{qo,qo,x) = ^—T^ Pi{x) , (2.13) 

where Ti(qo) = sin ^liQo)- The quantity Si{qo) is the phase-shift for a given angular 

momentum I and is determined in a standard manner. In Fig. 1 we show Re T{qQ, go, x, E) at 
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300 and 800 MeV laboratory energies together with partial wave sums up to a given angular 
momentum /. Note that E = Eiab/2. The strong peak of Re T{qQ,qQ,x, E) in forward 
direction requires high orders of Legendre polynomials for a correct description. This is of 
course especially pronounced for the higher energy. In Fig. 2 we display Im T{qQ,qQ,x, E) 
at the same energies together with its representation in partial wave sums. It can be seen 
that Im T{qQ, go, x, E) needs less partial wave amplitudes for its correct representation, the 
reason being that Im Ti is proportional to sin^ 5;, whereas Re Ti is proportional to cos 6i sin 5/. 
For large values of I the phase shifts become small, thus Im Ti decreases with 6f, whereas 
Re Ti only decreases proportional to 6i. 

An overview over the angular dependence of the full on-shell amplitude 
T{qo, go, x, E) as function of the energy is given in Fig. 3. Starting from a relatively fiat an- 
gular distribution at lower energies the peaking in forward direction develops with increasing 
energy. At the same time the angular range where the cross section is fiat and small increases 
with increasing energy, indicating that forward scattering dominates at higher energies. 

Next we consider the half-shell amplitude T{q,qQ,x, E) for two energies. Fig. 4 
shows that T(g, go, x, E) is rather small and structureless for all ofF-shell momenta g, with 
the exception of g being close to the on-shell momentum go. The most general amplitude, 
the fully off-shell amplitude T(g, g', x, E) is displayed in Figs. 5 and 6 (the real part is shown) 
as function of one off-shell momentum g and the angle x for two fixed momenta g'. Contrary 
to what one might expect, the strongest forward peaking does not occur for g being close to 
on-shell, but for g = q'. This agrees with the behavior of the driving term, which peaks for 
q = q' ■ We found this behavior for all energies E > 

All numerical and graphical examples considered so far refer to a potential of 
Malfiiet-Tjon type with repulsive and attractive parts (potential V'^^^ in Table I). Its strength 
is such that it supports a bound state at E =-2.23 MeV. Though this potential is of quite 
simple character, it captures essential features of NN interaction models based on meson 
exchange with respect to the propagator structure. 

III. THE OFF-SHELL T-MATRIX AT NEGATIVE ENERGIES 

In the context of Faddeev-Yakubovsky equations two-body t-matrices need to be 
evaluated at negative energies. In a three-body system the energy argument for the two-body 
t-matrix is given as E = Etot — -^(f 0- Here Etot is the total energy of the three particle 
system and -^q^ is the kinetic energy of the relative motion of the third particle with respect 
to the interacting pair, which is described by the t-matrix. For an interacting three-body 
system the relative momentum |q| is not conserved. Therefore it can have arbitrary values 
and E covers all energies below Etot- Thus we are interested to see, whether the angular 
dependence of the t-matrix evaluated at negative energies is similar to the one observed at 
positive energies. A second consideration is that bound states of the two-body system lead 
to poles in the t-matrix. The angular dependence at and around a pole should be dictated 
by the one of the bound state. We will investigate these questions and provide numerical 
illustrations. 

The formal solution to the Lippmann-Schwinger equation, Eq. ( p.l|) , is given by 

T{z) = V + V-^V, (3.1) 
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where H is the full two-body Hamiltonian. If this Hamiltonian supports a bound state 
at ^ = it follows immediately that 

T{z) v\<P,)^-r{h\V. (3.2) 

In momentum space representation Eq. ( p.2|) reads 

T(q',q,^) (q'|l/|0,)^_(0,|^|q). (3.3) 

The bound state obeys H\(j)b) = and has a certain fixed angular momentum 1, such 

that 

Since T is a scalar quantity, its behavior at and around the pole has to have the form 
T(q',q,^) '^'j:YiM')9i{Q')^^Yr^{q)g,{q) (3.5) 

m Z - J^b 

An z — hb 

^ RiWiQ,^ ■ q) 

z-Eb 

Here 

;>oo 

r7,(g)= / dq'q'\iq,q')(l)b,iiq') (3.6) 







2 /""^ 

vi{q,q') = -f drr'^ii{qr)V{r)ji{q'r). (3.7) 

TT JO 



with 



From Eq. (|3.6|) it can be clearly seen that the angular dependence of T exhibits 
a very characteristic behavior in the vicinity of the bound state poles, which is given by 
the Legendre function corresponding to the angular quantum number of the bound state. 
In order to illustrate the pole behavior, we choose the potential V^^^'^ of Table I, which 
supports a s-wave bound state at Es= -190.16 MeV and a p-wave bound state at Ep= 
-14.629 MeV. These binding energies are determined in a standard manner solving the 
Schrodinger equation for a fixed angular momentum. This is a simple, one-dimensional 
problem, whose solution also provides the function gi{q) and thus the residue i?/(g',g, q' • 
q) from Eq. (|3.6| ). The values of the binding energies can also be obtained by solving 
the two-dimensional integral equation Eq. (|2.8| ) and determining the pole position from 
the solution. Choosing the same integration points q in the partial wave projected, one- 
dimensional form and the two-dimensional form, the bound state energies Eb agree very 
well with the pole positions Epoi. For example, for 40 q-points (and 32 x-points) we find 
Eb{l = 0)= -190.162 MeV which has to be compared to Ep^ill = 0)= -190.164 MeV. Similarly, 
we find Eb{l = 1)= -14.6296 MeV compared to Epoi{l = 1)= -14.6296 MeV. These results can 
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be pushed to higher accuracy if desired. We also determine the residues at each pole from 
the solution of the two-dimensional integral equation and illustrate our result in Table II 
for the arbitrary choice of q' = q = qo = \Jm\E\ and the angle averaged quantity Ti = 
^ j\ dxPi{x)T{qQ, qo, x, E){E — Epoi{l)). As demonstrated in Table II, we approach the poles 
from both sides and the numbers closest to the poles agree very well with the corresponding 
residues calculated directly from the partial wave projected problem. 

The angular dependence of T{qo,qQ,x,E) with go = \f>T^]E\ for energies close 
to the two poles is displayed in Fig. 7. Both figures show that near and at the pole the 
t-matrix shows the characteristic behavior of the Legendre function associated with the 
angular momentum quantum number of the corresponding bound state. In Fig. 8 we show 
the angular dependence of T(go, qo, E) in the whole energy range around and in between 
the bound state poles. In order to include both poles we show {E — Es) {E — Ep)T {qo, qo, x, E) 
in Figs. 7 and 8. Starting from very small values of \E\, the angular dependence is first of 
s-wave character, then turns to a p-wave shape at and near the p-wave pole and then 
develops the forward peak known from corresponding positive energies. Note that due to 
the multiplicative factors, the peak turns upward in the two figures. When \E\ reaches the 
s-wave pole, T turns back to the pure s-wave behavior and then finally flips back into a 
strong forward peak. 

The fact that the angular dependence at negative energies is reminiscent of that 
at the corresponding positive energies, except for the characteristic behavior near the poles, 
leads us to suspect that the real parts of T might be quite similar to each other at energies of 
equal magnitude. This turns out to be the case as demonstrated in Fig. 9, where we display 
Re T[qo, qo, x, E) for the potential l^^^^ for different values \E\. It should be noted that the 
equality of Re T{qo, qo, x, E) for positive and negative energies is not the trivial consequence 
of ReT V, which does not hold. In order to demonstrate that V is significantly different 
from ReT, we also display V in Fig. 9. Comparing Re T{qo, qo, x, E) at the different energies, 
we have to conclude that the rescattering terms, ReT — V , behave more and more similar 
to each other for the same absolute values of the energy. 

Finally in Fig. 10 we display the real part of the fully off-shell t-matrix 
Re T{q, q' , x, E) as function of q and x for fixed energy E = 200 MeV and fixed momenta 
q' = 250 MeV/c and q' = 1000 MeV/c. As in Figs. 5 and 6 for positive energies Re T is 
most strongly peaked at g = q', which can be expected once the information of Figs. 5 and 
6 is known. 

IV. THE OFF-SHELL T-MATRIX IN THE SECOND ENERGY SHEET 

Two-body t-matrices might exhibit resonant behavior at positive energies or show 
a strong energy dependence near E = due to a virtual state. This latter case is realized 
for instance in the NN system for the partial wave state ^Sq. Our goal is to locate those 
resonances in the second energy sheet and investigate the characteristic angular dependence 
connected with the resonance or virtual state. 

The transition to the second energy sheet requires an analytic continuation of the 
Lippmann-Schwinger equation, Eq. ( ^2.21 ), into the second energy sheet, which we briefly 
describe here 0. For a complex energy z located in the upper half plane, we modify the 
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integration path as indicated in Fig. 11. The contribution along the closed path II gives a 
residue and we find 

T(q', q, z) = y(q', q) - nrmq, j dq"V{c{, g"g.)T(q"g„ q, z) 

ciVV(q',q")^-77rT(q",q,^). (4.1) 



Here = y/rnz, and the symbol at the second integral indicates the deformed integration 
path I. Since we deformed the integration path such that it is located above the energy z, 
we are able to take z into the lower half of the complex plane without hitting a singularity 
in the propagator. Once the energy z is located in the lower half plane, we can return with 
the integration path I to the real axis and have instead of Eq. (147 



T{q!, q, z) = 1/(q', q) - inmq^ J c/g'V(q', q"g^)T(q"g^, q, z) 

+ J d'q"Vi^, q")_J-^T(q", q, z). (4.2) 

m 

This equation is valid on the second energy sheet, which is reached from the upper rim 
of the cut along the positive real energy axis in the physical sheet. Due to the additional 
imaginary term, Eq. ( f4.2| ) has to be supplemented by another equation, which we obtain by 
choosing q' = qg^, 

T(q'g^, q, z) = V(q'g^, q) - iinnq, J dq"V{c(q^, ^'qz)T{q^'q^, q, z) 

+ J dYVic^q,, q")^^r(q", q, z). (4.3) 

m 

Similar to the bound state, which induces a nontrivial solution for the homogeneous 
equation related to Eq. (pTip, the homogeneous set of equations related to Eqs. ( |4.2| ) and 
( [4.3| ) has a nontrivial solution at discrete values of z. These discrete values either correspond 
to resonances with Re z > and Im z < or to virtual states with Re z < and 
Im z = 0. The fact that this homogeneous set of equations has a nontrivial solution 
together with the compactness property of the integral kernel means that T{z) has a pole 
at those energies z. We are interested not only in determining the positions of those poles 
but also in understanding the residues and their angular dependence. 

As it is obvious from Eq. ( [4.3| ), the driving term singles out the first entry as 
complex number. In order to simplify the formal steps leading to a determination of the 
residue, it is convenient to supplement the set of equations given in Eqs. ( [4.2| ) and ( |4.3| ) by 
another set in which the driving term has a complex entry in the second argument. The two 
sets can then be combined using the following matrix notation 

T(q',q,z) T(q',q'g„z) \ ^/ \/(q',q) ^q'^q'?.) 
,r(q'g^,q,2;) T{c(q,,qq^, z) ) VV(q'g^,q) V(q'g^,qg, 



(4.4) 



+ 



/rf3g'V(q',q") gf/rfg'V(q',q"g,) 
/ c?3g'T(q'g„ q") ql J dq"V{q'q,, q"g. 



T{q",q,z) T(q",qg„z) 
Tici"qz,q,z) T{q"q^,qq^,z) 
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Introducing the appropriate matrices, we write Eq. (4.4) as 

T = V + VGt. (4.5) 

We also need to study the corresponding homogeneous problem, which we want 
to write in the following form 

\{z)x = VGx- (4.6) 
In this form the eigenvalue is \{z) and the energy 2; is a parameter. Since VG is a compact 



operator, there is a discrete set of eigenvalues, which accumulate at \{z) = ||T0[- The 
physical resonances occur at those values z = Ej-es, for which X{Eres) = 1- 
In the following we choose z = E^es- Then we have 

X = VGx- (4.7) 
Since the kernel is nonsymmetric, we also have to consider the left hand eigenvalue problem 

= e^vG. (4.8) 

Defining 

$^ = 0^V (4.9) 

we deduce 

$^ = $^GV (4.10) 

or 

^ = V^G^. (4.11) 

Since V"^ = V, we obtain $ = 

In the immediate neighborhood of 2; = E^es and as a consequence of 

t(z) = (1 - VG)-iV (4.12) 

one has 

Tiz) ^(i_A(z))-il0^V, (4.13) 
where is a normalization factor. In the neighborhood of 2; = Eres we can put 

\{Z) ^ 1 + \'{Z) |,=E_ {Z - Eres) (4.14) 



and obtain 
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For the case of a bound state pole it is easy to prove that ^,^^.^| ^ ^ — jj: = 1 for a normahzed 

bound state \(f)b) and |x) = V\(j)b)- If we consider the right hand side of Eq. (|4.15| ) as function 
of an auxihary strength factor to the potential, we can adopt the normalization of the bound 
state and define 

f{z) X^^X^- (4.16) 

The final remark concerns the scalar nature of T(q', q, z). Since a resonant state 
has a unique angular momentum /, the function x(ct) will have the form 

X(q) = X/(g)>/m(q), (4.17) 

and we have to conclude that 

t(q',q,^) -^Xi{q') ^Xi{q)m-^). (4.18) 

If ReEp > and if ImEp < sufficiently small, the the t-matrix will feel the nearby pole 
also for real, positive energies z, and a resonance will occur in the differential cross section. 
In the case of a virtual state, like for ^Sq in NN scattering, the pole is located at ReEp < 
and ImEp = 0. For sufficiently small values of \ReEp\ the t-matrix will be strongly enhanced 
near and at z = 0. 

For our numerical realization we rewrite Eqs. ( f4.2|) and (^4.3| ) analogously to 
Eq. as 

T{q',q,x') = ^v{q',q,x, 1) 

POO rl 1 

+ / dq"q"' dx"v{q',q",x',x")^-^Tiq",q,x") 

Jo J -I ;z — 2— 

m 

-inmq, J ^ dx"v{q\ q^, x' , x")T{q^, q, x") (4.19) 



and 



T{q„q,x') = ^v{q„q,x,l) 



POO pi 1 

+ / dq"q"' dx"v{q,,q",x',x")^-j^Tiq",q,x") 
Jo J-i z — 2— 

m 

-iTimq^ j ^dx"v{qz,qz,x' ,x")T{qz,q,x"). (4.20) 



Here z = E = \E\e'^'^ with (j) < and = VrnE = y m|i?|e*'^/^. The nontrivial solution 
to the homogeneous system of Eq. (|4.7|) has a fixed angular momentum. When employing 
Eq. ( [4.17| ), we obtain 

POO 2_ 

Xi{q) = dq'q'^vi{q,q') -J2Xi{q') - iT^mq,Vi{q, q,)xi{qz) (4.21) 

Jo z - ^ 
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and 



dq'q' vi{q^,q') 



-xiiq') - iT^mq^viiqz, qz)xiiqz 



(4.22) 



We used the above given equations Eqs. ( [4.21[ ) and ( |4.22| ) to determine the location of the 
resonances in the second energy sheet. For s- waves they are usually called virtual states and 
are located on the negative energy axis. For partial waves with / = 1 or higher, the energy 
eigenvalues have a positive real part and a negative imaginary part. For varying potential 
strength, they move along trajectories in the complex energy plane. In the following, we 
numerically study two different cases, namely a s-wave virtual state supported by potential 
model V^'^'^ of Table I and a p-wave resonance of potential model V^^^^ of Table I. The 
corresponding trajectories are listed in Table III and shown in Fig. 12. 

For the case of the s-wave virtual state we started from potential model V^^^ of 
Table I. The s-wave phase shift of this potential has an effective range expansion with a 
scattering length = —23.5818 fm and an effective range = 2.8789 fm. Using these 
values, the pole position for the S-matrix can be estimated via the effective range expansion 
as 



q: 









1 


i 




/ f 








V rs\as 


^ s . 



(4.23) 



which leads in our specific case to a position of the virtual state E^^^ = ^ 
—0.06738 MeV, where we used m = 938.9 MeV. This number is very close to the exactly 
calculated value given as Ep{l = 0) = —0.06663 MeV, which is given in Table III. In the 
same table we quantify the s-wave trajectory as a function of the strength parameter Va of 
potential V^^\ 

For obtaining a p-wave resonance state, we start from potential V^^^^^ of Table I, which 
supports a p-wave bound state, and decrease the attraction by decreasing the strength 
parameter Va until the bound state turns into a resonance state. Selected values for the so 
obtained trajectory for the p-wave resonance are listed in Table III. Of course, this model 
does not correspond to the reality of an NN force, even for the lowest value of = 4.7 
given in Table III, the binding energy of the s-wave bound state is still Eg = —52.52 MeV. 
Nevertheless, this example illustrates in a simple manner, what can be expected for other 
cases like an effective nucleon-nucleus interaction, which supports low energy resonances for 
certain angular momentum states. Qualitatively the same picture would emerge. 

We would like to mention that we solved the homogeneous set of Eqs. ( |4.21| ) and 
( |4.22| ) by the very efficient power method Regarding the notation of Eq. ([4.7|) , one 

has to determine the eigenvalue X{z) and vary the energy such that X{z) = 1. For the 
potentials used here, the largest eigenvalue in magnitude was always an unphysical one with 
a negative real part generated by the repulsive short range piece of the force. Once the 
largest eigenvalue is determined, we introduce a new integral kernel, consisting of the old 
one minus that specific eigenvalue. The in this way defined new kernel has then the physical 
eigenvalue as the largest one in magnitude. 

Next we investigate the solution of the t-matrix in the second energy sheet as 
given by Eqs. ( [4.19| ) and (4.20). We are interested in verifying the location of the pole 
as well as the angular dependence of the residue at the pole as given in Eq. ( [4.18|) . We 
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illustrate our findings for the potential V^^'^^ of Table I, which has a p-wave resonant state 
at Eres{l = 1) = (4.3177 - 22.2386) MeV and for the potential V^^'^ of Table I, which has a 
virtual s-wave state at Eres{l = 0) = —0.06663 MeV. For the p-wave resonance we show in 
Figs. 13 and 14 the angular dependence of Re [{E — Eresi} = ^))T{qE,(lo,x, E)] as function 
of the complex energy E located along two straight lines going through the pole position. 
In Fig. 13 we start on the real axis at = 4.3 MeV and successively increase the imaginary 
part of E. For E = (4.3 — z2.4) MeV we clearly see an angular dependence characteristic of 
a p-wave residue. Since the width of the resonance is relatively small, the p-wave behavior 
is present along the whole vertical energy line including the point on the real axis. For 
our second choice of energy line (Fig. 14), a declined line starting from zero energy, starts 
with a behavior being a mixture of s- and p-wave, but still relatively fiat. Approaching the 
resonance, the shape becomes predominantly the one of a Pi{x) given by the p-wave residue. 

A corresponding study based on the potential model V^^'^ is shown in Fig. 15 
for the negative real axis in the second sheet, where the neighborhood of the virtual pole 
position is considered. In the vicinity of the virtual state, the residue exhibits perfect s-wave 
characteristics. 

Finally, we would like to demonstrate the structure of the t-matrix as given 
in Eq. ( [4.18| ) in a numerical example. First, we numerically verify that lim^^^;,,^^ (z — 
Eres)T{q',q,x,z) behaves like Pi{z). Then we determine Xiil) by comparing the numer- 
ically calculated quantity {z — Eres)T{q', q, x, z) to the form given in Eq. ( |4.18| ). Instead of 
dividing by Pi(x) we use 



Xi{q)= lim i/27r / dxPi{x){z ~ Eres)T{q,q,x, z). {AM) 

The values of xi (q) for z approaching Eres are shown for a few arbitrarily selected momentum 
points q in Table IV. They stabilize for z approaching the pole position Eres = (4.318 — 
z2.239) MeV. The last row in Table IV shows the values obtained from the solution of the 
homogeneous set of Eqs. ( |4.21| ) and ( |4.22| ). They have been normalized at one q point to 
the Xi{l) extracted from T{q,q,x,z Eres)- The agreement is perfect. 

Finally, we directly verified the separable structure of T(g', q,x,z ^ Eres) as given 
in Eq. (|4.18|) by evaluating T = 2it J\ dxPi{x){E — Eres)T{q' , q, x, E) very close to ii^ = Eres 
for different q' ^ q and comparing them to the values obtained via Table IV. The agreement 
is again perfect and the values for a selected set of momentum points q' and q are given in 
Table V. 

V. SUMMARY 

Two nucleon scattering at intermediate energies of a few hundred MeV requires 
quite a few angular momentum states in order to achieve convergence of e.g. scattering 
observables. This is even more true for the scattering of three or more nucleons upon 
each other. An alternative approach to the conventional one, which is based on angular 
momentum decomposition, is to work directly with momentum vectors, specifically with the 
magnitudes of the momenta and the angles between them. We formulated and numerically 
illustrated this alternative approach for the case of two-body scattering, which includes the 
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approach towards bound states and resonances. The angular dependence of the two-body 
t-matrix is directly determined from the Lippmann-Schwinger equation, which now is a two- 
dimensional integral equation in contrast to the one-dimensional one for a fixed angular 
momentum in a partial wave formulation. This two-dimensional integral equation is quite 
easily numerically tractable. We determined the angular dependence of the on-shell, half- 
shell and fully off-shell t-matrix as function of the scattering energy and different choices 
of momenta. As two-body force we concentrated on a superposition of an attractive and 
repulsive Yukawa interaction, which is typical for nuclear physics. We neglected spin degrees 
of freedom in all our studies. 

We want to briefly summarize our results. The on-shell t-matrix develops a strong 
forward peak as the energy increases, which is more and more difficult to build up in a 
calculation based on angular momentum decomposition, but relatively simple accessible in 
our approach using momentum vectors. The angular dependence of a half-shell t-matrix 
is strong only around the on-shell momentum and rather mild otherwise. For a fully off- 
shell t-matrix T{q, q', x, E) the strong angular dependence occurs for q — q', which do not 
necessarily have to coincide with the on-shell momentum. At negative energies the t-matrix 
has poles located at the bound state energies, if those exist. As example we investigated 
s- and p-wave bound states. The numerically determined t-matrix turned out to be very 
well under control even quite close to the the bound state poles, where the homogeneous 
version of the Lippmann-Schwinger equation has a nontrivial solution. We determined the 
angular dependence of the t-matrix at the two poles, at energies between them and at 
energies way below the deepest bound state. Directly at the poles the angular behavior 
displays the characteristics of the Legendre polynomial of the same angular momentum / 
as the bound state. Between the poles as well as for energies below the last bound state 
the t-matrix exhibits the same forward peaking as visible at positive energies. This latter 
result is interesting by itself. The angular dependence at positive and negative energies is 
very similar. More quantitatively, we found that the real parts of the t-matrix are extremely 
close to each other at positive and negative energies of equal magnitude long before this 
statement becomes trivial due to the validity of the Born approximation T = V. 

Finally we studied the analytical continuation of the Lippmann-Schwinger equation 
into the second energy sheet, which is reached through the cut along the positive real axis 
of the physical sheet. In the lower half plane possible resonance poles of the t-matrix are 
located, which are of course of interest only if they are close to the real axis. As example 
we studied a p-wave resonance and mapped out its pole trajectory by varying the potential 
strength. At the pole, the t-matrix achieves a separable form, which we verified numerically. 
We also verified the characteristic angular dependence of the t-matrix close to the resonance, 
which corresponds to the Legendre polynomial of degree / of the angular momentum state of 
the resonance. For negative energies in the second sheet we investigated the pole trajectory 
of a virtual s-wave state, which is of interest in the NN system for the quantum number ^Sq. 

Summarizing we can state that the two-dimensional Lippmann-Schwinger equa- 
tion can be handled quite easily in a numerically very reliable manner. In this approach 
one determines directly the angular dependence of the t-matrix for arbitrary momenta and 
energies. Once supplemented by spin degrees of freedom this approach will be of interest 
in the NN system. In addition, this approach will be generalizable to systems with three 
particles, like three nucleons or two nucleons and a meson. In the case of three nucleons. 
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Faddeev calculations at e.g. 150 MeV and higher are getting quite tedious because of the 
very many orbital angular momentum states involved ^ and a direct, three-dimensional 
approach appears to be preferable. First steps in this direction are under way. 
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TABLES 



TABLE I. Parameters of the Malfliet-Tjon type potentials. As conversion factor we use units 



such that nc=197.3286 MeVfm=l. 





Va 


MA [MeV] 


Vr 


m [MeV] 




3.1769 


305.86 


7.291 


613.69 








10.0 


2000. 




6.0 


305.86 


7.291 


613.69 




5.1 


305.86 


7.291 


613.69 




2.6047 


305.86 


7.291 


613.69 



TABLE II. Determination of 7] = 1/q dxPi(x)T{qQ, qq, x, E){E — Epoi{l)) as function of E 
close to the I = 1 and / = poles. The values for the constants are co=2 and ci=2/3. The entry 
p.w. indicates the value determined from the partial wave projected problem. 



E [McV] 


fi{E) [MeV fm2] 


E [MeV] 


fo{E) [MeV fm2] 


-14.60 


1.08424 


-189.8 


4.37911 


-14.61 


1.08513 


-189.9 


4.37847 


-14.62 


1.08597 


-190.0 


4.37772 


-14.63 


1.08588 


-190.1 


4.37722 


-14.64 


1.08778 


-190.2 


4.37693 


-14.65 


1.08857 


-190.3 


4.37608 


p.w. 


1.08684 


p.w. 


4.37685 



TABLE III. Pole trajectories in the complex energy plane as function of the strength parameter 
Va for the potentials V^^^^^ (p-wave pole trajectory) and y(^) (s-wave virtual state trajectory). 



Va {V('^'^) 


Eresil = 1) [MeV] 


Va (F(^')) 


Eresil = 0) [MeV] 


5.4 


0.9267 -i 0.1821 


2.6047 


-0.06663 


5.3 


2.1972 -i 0.7101 


2.6 


-0.07380 


5.2 


3.3254 -i 1.4134 


2.5 


-0.31004 


5.1 


4.3177 -i 2.2386 


2.4 


-0.69870 


5.0 


5.1801 -i 3.1552 


2.3 


-1.22970 


4.9 


5.9178 -i 4.1413 


2.2 


-1.89209 


4.8 


6.5358 -i 5.1806 


2.1 


-2.67878 


4.7 


7.0388 -i 6.2597 


2.0 


-3.57733 



16 



TABLE IV. Determination of xi{q) ™ the vicinity of the p-wave resonance of potential model 
V^^'^^ as function of the complex energy E. 



E [MeV] 


Xi(q=94.89 MeV/c) 


Xi(q=292.73 MeV/c) xi(q=527.50 MeV/c) 


;j^i(q=758.09 MeV/c 


4.32 -i 2.0 


1.105 -i 0.400 


1.856 -i 0.729 


1.149 -i 0.473 


0.491 -i 0.212 


4.32 -i 2.10 


1.103 -i 0.400 


1.851 -i 0.729 


1.144 -i 0.471 


0.488 -i 0.207 


4.32 -i 2.20 


1.101 -i 0.400 


1.846 -i 0.729 


1.139 -i 0.468 


0.484 -i 0.202 


4.32 -i 2.23 


1.101 -i 0.400 


1.844 -i 0.729 


1.138 -i 0.467 


0.483 -i 0.201 


4.32 -i 2.4 


1.098 -i 0.400 


1.835 -i 0.729 


1.129 -i 0.461 


0.477 -i 0.192 


4.32 -i 2.5 


1.096 -i 0.400 


1.830 -i 0.729 


1.124 -i 0.459 


0.473 -i 0.187 


p.w. 


1.102 -i 0.400 


1.844 -i 0.729 


1.138 -i 0.468 


0.482 -i 0.200 



TABLE V. Comparison of T = 27t j\dxPi{x){E - Eres)T{q,q',x,E) with Xi(g)xi(g') for a 
fixed value q' and different values of q at the p-wavc resonance of potential model The 
energy E for calculating f was fixed at E = (4.32 - i2.23) MeV. 



q [McV/c] 


q' [MeV/c] 


T[MeVfm^] 


xMximMeVfm^] 


292.73 


94.89 


1.7398 -i 1.5395 


1.7397 -i 1.5398 


527.50 


94.89 


1.0668 -i 0.9683 


1.0667 -i 0.9689 


758.09 


94.89 


0.4522 -i 0.4134 


0.4521 -i 0.4140 
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FIGURES 



FIG. 1. Angular dependence for the real part of the on-shell t-matrix, Re T{qo,qo,x,E). At 
Eiab = 300 MeV (a) the dashed line represents the partial wave sum up to j = 2, the dash-dotted 
line the sum up to j = 4 and the solid line the sum up to j = 6. The solid bullets stand for the 3D 
calculation. At Eiab = 800 MeV (b) the dashed line represents the partial wave sum up to j = 6, 
the dash-dotted line stands for the sum up to j = 9, and the solid line for the sum up to j = 12. 
Again, the solid bullets represent the 3D calculation. 

FIG. 2. Angular dependence for the imaginary part of the on-shell t-matrix, ImT{qo, qo, x, E). 
At Eiab = 300 MeV (a) the dashed line represents the partial wave sum up to j = 2, the solid line 
the sum up to j = 4. The solid bullets stand for the 3D calculation. At Eiab = 800 MeV (b) the 
dashed line represents the partial wave sum up to j = 3, the solid line the sum up to j = 6. The 
bullets stand for the 3D calculation. 

FIG. 3. The angular dependence of the on-shell t-matrix as function of the laboratory energy 
from Eiab = 50 MeV to Eiab = 1000 MeV. On the left ReT{qo, qo, x, Eiab) is displayed, on the right 
ImT{qo,qo,x,Eiab). 

FIG. 4. The angular dependence for the real part of the half-shell t-matrices. Re T{q, qo, x, E), 
is displayed for Eiab = 200 MeV (a) and Eiab = 500 MeV (b). 

FIG. 5. The angular dependence for the real part of the off-shell t-matrix, Re 
T{q,q' = 2hQMeV/c,x,E), is displayed for E^b = 400 MeV. 

FIG. 6. The angular dependence for the real part of the off-shell t-matrix, Re 
T{q,q' = 1000Mey/c,x,£;), is displayed for Eiab = 400 MeV. 

FIG. 7. The angular dependence for {E — Es){E — Ep)T{qo, qo, x, E) for energies E around the 
p-wave pole (a) and the s-wave pole (b). 

FIG. 8. The angular dependence for {E — Es){E — Ep)T{qo, qo,x, E) as function of the energy 
from E=-200 MeV to E=-l MeV. Note the characteristic angular behavior around the p- and 
s-wave poles as well as the strong forward peak between the two poles and below the s-wave pole. 

FIG. 9. The angular dependence for the real part of the on-shell t-matrix Re T{qo, qo, x, E) for 
|£^|=200, 400 and 800 MeV. For comparison, the angular dependence of the driving term, , is 
also shown as dash-dotted line. 

FIG. 10. The angular dependence for the real part of the off-shell t-matrix Re 
T{q,q',x,E = 2mMeV) is displayed for q'=250 MeV/c (a) and for q'=1000 MeV/c (b). 
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FIG. 11. Modified integration path in the complex q plane for the analytic continuation of the 
Lippmann-Schwinger equation into the second energy sheet as described in the text. 

FIG. 12. The pole trajectories of the p-wave and s-wave pole in the second energy sheet as 
given by the potentials and F^^) of Table I. 

FIG. 13. The angular dependence for the real part of {E — Eres)T (qs , Qo, x, E) as function of the 
complex energy E around the p-wave resonance of potential model y (^^) . The off-shell momentum 
go is fixed at 100 MeV/c. 

FIG. 14. Same as Fig. 13, but for a different path of the complex energy E. 

FIG. 15. The angular dependence of {E — Eres)T{qE, QO: x, E) as function of the negative energy 
E in the second energy sheet around the s-wave virtual state of potential model V^^\ The off-shell 
momentum qo is fixed at 100 MeV/c. 
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